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TOWARDS UNDERSTANDING THE MECHANISM OF RECEPTIVITY AND BYPASS 
DYNAMICS IN LAMINAR BOUNDARY LAYERS* 


D.G. LASSEIGNE+, W.O. CRIMINALS*, R.D. JOSLIN§, AND T.L. JACKSON* 

Abstract. Three problems concerning laminar-turbulent transition are addressed by solving a series of 
initial value problems. The first problem is the calculation of resonance within the continuous spectrum of 
the Blasius boundary layer. The second is calculation of the growth of Tollmien-Schlichting waves that are 
a direct result of disturbances that only lie outside of the boundary layer. And, the third problem is the 
calculation of non-parallel effects. Together, these problems represent a unified approach to the study of 
freestream disturbance effects that could lead to transition. Solutions to the temporal, initial-value problem 
with an inhomogeneous forcing term imposed upon the flow is sought. By solving a series of problems, it is 
shown that: 

• A transient disturbance lying completely outside of the boundary layer (‘an lead to the growth of an 
unstable Tollmien-Schlichting wave. 

• A resonance with the continuous spectrum leads to strong amplification that may provide a mecha- 
nism for bypass transition once nonlinear effects are considered. 

• A disturbance with a very weak unstable Tollmien-Schlichting wave can lead to a much stronger 
Tollmien-Schlichting wave downstream, if the original disturbance has a significant portion of its 
energy in the continuum inodes. 

Key words, boundary layer, receptivity 

Subject classification. Fluid Dynamics 

1. Introduction. In previous work [3], [8], the authors have shown a strong correlation to the solution 
of a temporal, three-dimensional, initial-value problem and the direct numerical simulation of the spatial 
problem. The methodology consisted of solving the linear disturbance equations subject to a series of initial 
values. These solutions are relatively easy, fast, and inexpensive to calculate. The corresponding spatially 
evolving flow was then determined by direct numerical simulation using the full Navier- Stokes equations 
and the two solutions were compared. During the period of transient growth for both channel flow and the 
laminar boundary layer, the two approaches agree quite well. Thus, it is reasonable to use the inexpensive 
and fast solutions of the temporal, initial-value problem as a means to conduct numerical experiments that 
can lead to greater understanding of receptivity and bypass mechanisms. Of course, the suitability of this 
approach must be continually confirmed by selectively using the more expensive direct numerical simulation 
to compare with the major results. For channel flow [3], the dynamics of specific initial conditions were 
determined, and the growth of disturbance energy compared to initial conditions which produce optimal 
growth of disturbance energy. The optimal initial conditions were determined by appropriately expanding 
the initial condition, finding the solution to a relatively few number of initial value problems, and determining 
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the coefficients which maximized the disturbance energy. The transient period of the laminar boundary 
layer [8] was investigated in a similar manner. By using this approach, the contributions of the continuum 
modes are properly assessed, and it is determined that they must play a crucial role in the analysis of 
bypass mechanisms. It is only natural that we extend our studies to investigate the effects that freestream 
disturbances have on the laminar boundary layer. 

Receptivity has traditionally been divided into two broad categories, forced receptivity or natural re- 
ceptivity, based on the physical and mathematical descriptions. Forced receptivity is characterized by the 
experiments of Nishioka and Morkovin [9] where disturbances of limited spatial extent are introduced in 
the freestream downstream of the leading edge. The case of natural receptivity is characterized by the 
experiments of Boiko, Le.,[ 1] where a disturbance field upstream of the leading edge of a smooth plate is 
generated. Great care is used to establish a zero pressure gradient at the plate, so that the mean flow (in 
the absence of the grid) has a Blasius profile. By assuming that disturbances are kept at a level in which 
linear theory applies, the differences between the two types of receptivity are seen mathematically in that 
the forced receptivity problem is governed by a set of inhomogeneous linear partial differential equations in 
time and space and natural receptivity by the homogeneous problem. Downstream of the imposed distur- 
bance, the forced receptivity problem is the same mathematically as the natural receptivity problem. In 
addition to the two traditional categories of receptivity, there is a third case of receptivity that results in 
either an inhomogeneous problem or a homogeneous problem depending on the modeling approach taken. 
This case concerns the scattering of freestream disturbances by localized surface irregularities (e.g. acoustic 
disturbance or free stream turbulence with surface roughness or surface blowing and suction). Viewed as a 
perturbational problem, an inhomogeneous problem results; however, viewed as a problem of changing mean 
How, a homogeneous problem results. Here, it is referred to as the naturally-forced receptivity problem since 
it has characteristics of both of the traditional categories. 

The procedure used here, integrating the linear disturbance equations of temporal stability theory as an 
initial value problem, is straightforward and simple. It has already been demonstrated to be able to agree 
well with tin 1 direct numerical simulation of the spatial problem since in every numerical calculation, the 
complete solution including the continuum eigenfunctions of the Orr-Sommerfeld problem and all discrete 
modes is determined. Only afterwards is this solution interpreted in terms of the individual modes of the 
Orr-Sommerfeld equations. The theory of expanding the solution in terms of Orr-Sommerfeld eigenfunctions 
is presented in the classic work of Salwen and Grosch [10] while the description of the important continuum 
modes is found in Grosch and Salwen [6]. 

The rest of the paper is organized in the following manner. In Section 2, the equations governing the 
evolution of a disturbance under the assumptions of parallel linear theory and subjected to forcing terms 
are presented. The forcing terms can be interpreted as resulting from the forced receptivity problem or the 
nonlinear interaction of two linearly independent disturbances in the naturally forced receptivity problem. 
Fourier transforms in the streamwise and spanwise direction reduce the equations to partial differential 
equations in time and the vertical variable y. In the case of naturally forced receptivity, it is determined that 
the only part of the force that can affect the normal velocity component, and therefore the generation of 
Tollmien-Schlichting waves, is the divergence free component of the force. In Section 3, two major issues are 
explored by constructing the complete solutions of the temporal inhomogeneous problem using a model mean 
flow. The first issue deals with the difference between specifying the forcing by prescribing the individual 
forcing components in the momentum equations (naturally forced receptivity) or by specifying a vorticity 
source (forced receptivity). It is shown that by specifying the vorticity source, there is an immediate response 
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throughout the boundary layer which leads to a much higher degree of receptivity than the specification of 
forcing in the momentum equations. The second issue deals with the possibility of resonant forcing. Both 
these important issues are later addressed by numerical calculations using the continuous mean value profile, 
and the analytical results of this section are verified. In Section 4, the construction of the solution to the 
inhomogeneous problem in terms of all of the eigenfunctions of the Orr-Sornmerfeld problem is presented. 
The gain in the coefficients of an eigenfunction expansion of the solution from before the imposition of 
forcing to after the imposition of the forcing is used to characterize the receptivity due to the forcing. Most 
importantly though, it is noted that the results after the imposition of the disturbance should be considered 
as initial values at a higher Reynolds number. Thus the transfer of the solution in terms of an Orr-Sommerfeld 
eigenfunction from one Reynolds number to another is presented. It is also shown in this section, that if 
the forcing is specified as forcing components in the momentum equations, that the eigenfunction expansion 
of the solution to the inhomogeneous problem has a simple form. From this solution, it is seen that the 
resonance found in Section 3 is indeed a resonance of the forcing with the continuum modes. All results 
presented in Section 4 are verified by numerical calculations. 

2. Basic Governing Equations. For the fiat-plate boundary layer, the fluid is taken as oik 1 of constant 
density with the basic flow approximated as parallel with U = U(y), V — IF = 0. The instantaneous How is 
decomposed into a basic state, (U, V , \\\ P), plus a time-dependent disturbance to this basic state, (u,v,w,p). 
Then, the nondimensional linearized equations of motion are written as 


(2.1) 


du dv dw 
dx dy dz ' 


(2.2) 

(2.3) 

and 

(2.4) 


du du dU dp , 

m +l ^ + d^ v + ^ = R 


d 2 u d 2 u d 2 u 
dx 2 + dy 2 + dz 2 


+ A, 


dt dx dy 


d 1 v d~v d 2 v 

d7 2 + dtf + di? 


+ D. 


dw rr dw dp „ 

li + lj oi*oi = R ^ 


d 2 w d 2 w d 2 w 

+ Th/ 2 + 


+ c. 


The length scale chosen for non-dimensionalization is the displacement thickness, and thus the non-dimensional 
mean velocity is related to the Blasius boundary layer solution by U (y) = f^^y) with /(?;) being the solution 
to 


(2.5) /'" + !//" = 0 subject to /(0) = /'(0) = 0, /'( oo) = 1. 

where <5,* = / 0 °° (1 - /' (tj)) dr] « 1.7208. The nondimensional mean profile satisfies (1 — U{y)) dy = 1 
and U (2.856) = .99 gives the outer edge of the boundary layer. Here, R is the conventional Reynolds number 
based on the displacement thickness, and time is nondimensionalized by the advective time scale. In addition, 
A, B and C are the nondimensional components of a prescribed force, F(x, y, z, t). 

Then, by using the Fourier transformations defined with respect to x and 2 as 


(2.6) 


v(a,y,&,t) 



v(x , y : z , t)e l{ax + fiz ^dxdz 
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(etc. for u,w, p, ,4, B and C), equations (1) to (4) become 


(2.7] 


-i(au + Pun 4- 7T- — 0, 


( 2 . 8 ) 




iaUu 4- 17 7 e — mp — R 


-l 


<9 2 ix 

tty 2 


— 7 2 u 


4- -4, 


(2.9) 
and 

(2.10) 


d v . dp _ j 

_- tabv + -_ R 


d' 2 v ' 
^ 


dw 

dt 


— iotUw — ipp = /? 


- j 


3 2 tc 


% 


- - 7 2j " 


+ B. 


4- C , 


respectively, where U' = dU/(ly and y 2 - a 2 4- /i 2 . The Squire transformation, written as 


(2.11) aw 4- flw — 7 u , 

(2.12) - -f me = 71&, 

and combined with operations on (2.7) to (2.10) enables the elimination of the pressure to obtain the pair 
of equations 


(2-13) 




An + iall n v — R 1 A An 4- i 


( dA dC 
dy dy 


-TB 


and 

(2.14) 


' d_ 
dt 


— iaU 


w = sin <t>U'v 4- R 1 Ate — sin 4>A + cos <pC , 


where A is the linear operator 
(2.15) 


A = 


d 2 

% 2 


Here, sin 0 = /J/7, At) is proportional to the difference of the vorticity components in the x — z plane 
( Ac = ipu)y - ?*qw : ) and w is proportional to the normal vorticity component (u) y = ijw). The first 
inhomogeneous term on the right hand side of (2.14) has been denoted as vortex tilting and acts as forcing 
for the normal vorticity. Such tilting is a product of the mean vorticity in the spanwise direction (Q z — -[/') 
and the perturbation strain rate (dv/dz). For a three dimensional disturbance, the vortex tilting gives 
rise to the increase of the normal vorticity. It is clear that the solutions of (2.13) and (2.14), combined 
with continuity and the Squire transformation, are equivalent to solving (2.7) to (2.10). Lastly, p can be 
determined once the v component of velocity is known. In either case, solutions of the equations are subject 
to imposed initial conditions and the following boundary conditions at the the plate, namely 


(2.16) 


dv 

v{0J) - —{0 J) = w(0,f) = 0, 
dy 
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as well as boundedness conditions in the freest ream. 

To evaluate the other velocity components, the quantities v and w are computed from (2.13) and (2.14), 
respectively. Them the Squire transformation, (2.11) and (2.12), is inverted to give 


(2.17) 
and 

(2.18) 


„ i cos 0 dv . _ 

u — — sin <p w 

7 OV 


„ i sin 6 dv 
w — — — h COS 0 w . 

7 9y 


An energy density in the (7, 0) plane as a function of time is defined as 


(2.19) 


E'(t) 



+ |j ;| 2 + |w| 2 ] dy 


or, in terms of the variables v and t», this becomes 


nOO 

(2.20) E*(t) = / 

Jo 

In using these definitions, the class of disturbances is more restricted than those that are only bounded at 
infinity. However, in all calculations that follow, both the specified forcing and the initial conditions are 
chosen as decaying as y -» 00 and therefore the energy integral defined as (2.20) converges for all time. The 
continuum is still represented in these calculations, as it is the integral over all the continuum modes, not 
an individual continuum mode that appears in the complete solution. The total energy of the perturbation 
can be found by integrating (2.20) over all 7 and 0. A normalized energy density, namely 


+ + ■ 


dv 

2 

dy 



dy . 


(2.21) £ «> = f£iV 

measures the growth in energy at time t for a prescribed initial condition at , t = 0 and subject to prescribed 
forcing. 

2.1. Interpretation of the Forcing: F. By knowing the perturbation velocity components, the 

perturbation vorticity is determined by appealing to their definitions, namely 


( 2 . 22 ) 




dw 

= -7T- + 10V, 

dy 


(2.23) 


LUy = —iflii + iaw = i/yw, 


and 

da 

(2.24) a;- = —iav — — , 

dy 

in wave space. The governing equations, (2.7) to (2.10), are recast into equations governing vorticity and 
pressure so that the effect the chosen forcing has on the flow can be interpreted. These equations are: 


(2.25) 


wr iM 


— iaU'w = R 


dC 

1 Au3j. + — — h i8B , 
oy 
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(2.26) 


I - iaV 


LJy - iSU'v 


R 1 Au>j, — if3A + iaC, 


(2.27) 


and 

(2.28) 


<9 TT 

di~ ,aL 


- U" v = R- 1 Ai) ; 


. - 

u*£, 

oy 


Ap = 2iaU'v — i(aA + ^C) + 




The components of the force F = (.4, £,C) are specified according to the particular problem at hand. 
Since any vector field can be decomposed into a field that is divergence free and a field that is curl free, 
and since the forcing terms in the vorticity and pressure equations are the curl and divergence of the force, 
respectively, this is a proper place to start. First, suppose that F represents the component of the force that 
is curl free, Le., 


(2.29) 


VxF = 0. 


In Fourier space this means 

dC * - - dA 

(2.30) — + i(3B = 0; aC - HA = 0; iaB + - 0. 

dy oy 

For this choice, the forcing terms appearing in the vorticity equations vanish, and only the pressure dis- 
turbance is directly affected. To specify this case, any one of the components of the force is considered as 
arbitrary while the other two components must satisfy (2.30). 

Next, suppose that F represents the component of the force that is divergence 4 free, z.e., 

(2.31) V ■ F = 0, 


or, in Fourier space, 

OB 

(2.32) -i(aA + &C) A = 0. 

Now, non-zero forcing terms appear only in the vorticity equations. To specify this case, two components of 
the force are considered as arbitrary while the last must satisfy (2.32). For convenience, the two independent 
modes of forcing in the normal velocity equation (2.13) without forcing in the normal vorticity equation 
(2.14) given by 

(2.33) S/0, and — {3 A + aC = 0 

and the opposite, forcing in the normal vorticity equation (2.14) without forcing in the normal velocity 
equation (2.13) given by 

(2.34) B - 0, and - 0A + aC = 

will be the only pair of divergence free forces considered. In the first instance, the forcing in (2.13) is the 
Laplacian of B. 
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3. Complete Solution to an Inhomogeneous Model Problem. It has been shown by Crimi- 
nale and Drazin [2] that the viscous boundary layer problem can be solved analytically for large values of 
the Reynolds number by modeling the mean flow. This technique is particularly appropriate for studying 
boundary-layer free stream response to external forcing. In this limit, an outer layer, i.e., the free stream, 
exists that approximates an irrotational flow of an inviscid fluid, an intermediate* layer exists that approxi- 
mates a rotational flow of an inviscid fluid, and a viscous sublayer exists at the plate. To leading order, the 
viscous terms are zero in the first two regions. The basic mean boundary layer flow is modeled as a uniform 
parallel stream above a uniform shear flow with a solid boundary below. Thus, for this entire section, the 
basic velocity is taken as 


(3-1) 


U(y) = (U, 0,0) 


r (i,o,o) 


{ (y, o,o) 


y > 1 


0 < y < 1. 


After substituting U from (3.1) and neglecting the viscous terms, the governing equations, (2.13) and (2.14), 
reduce to 


(3.2) 

and 

(3.3) 


' d_ 
dt 


iaU 


. dA DC 
Av = i «— + ji— 

V °V vy 


fB 


lot 


iaU 


iv = sin 0 U'v — sin 0 A T cos 0 C, 


which must be solved on either side of the interface located at y — 1 and coupled to the viscous sublayer 
below. Both the normal velocity v and the perturbation pressure 

(3.4) -j 2 p = 


j,- inU 


dv 

Qy 


+ iaU'v - ?’(ad -f flC) 


are required to be continuous at y — 1. The plate conditions where v = v* = w = 0 at y — 0 are secured by 
the solutions influenced by viscosity. 

In what follows, a number of issues concerning the specification of the forcing in the vorticity equation 
will be resolved, and results of using models based on these concepts agree with the numerical results using 
the Blasius profile. These issues are directly related to whether the inhomogeneous terms represent the forced 
receptivity problems or the nonlinear interaction terms in the naturally forced receptivity problem. In the 
latter case, the inhomogeneous terms should be consistent with the velocity field since they have evolved 
along with the flow. In the former case, the externally applied forcing does not have to be consistent with 
a velocity field and this results in generation of vorticity at the plate. First, however, the issue' of forcing in 
the pressure equation is examined. 


3.1. Forcing in the Pressure Equation. In considering forcing due to pressure, the component of 
a force, F, with curl equal to zero is all that is required. Therefore, for this analysis, it. is assumed that 
V x F = 0. The governing equations are 


(3.5) 
and 

(3.6) 


l~ iaU 


Av = 0. 


f s ,, 

si~ u 


w = sin 0 U'v 
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In like fashion, the equations for the vorticitv ((2.25)-(2.27)) have no forcing terms. In short, unless vorticity 
is an input initially, it cannot be generated and solutions to (3.5) and (3.6) have only irrotational portions. 
The pressure equation does, on the other hand, have an extra term due to forcing and is 

dB 

(3.7) Ap = 2 iaU'v - i(aA -f fiC) + - — . 

dy 

This relation can be cast in real space as well, namely 

(3.8) V 2 p = -2U'—-VF. 

Since V x Vp = 0 and defining 

(3.9) VP = Vp+F, 


(equation (3.8) becomes 

(3.10) V 2 P = _ 2t r/ — . 

Thus. P plays the same role as p when no external forcing in pressure is present. Seeing that this component of 
the forcing cannot generate vorticity, and therefore Tolhnien-Schlicting waves, at this order of appoximation, 
attention is now focused on the divergence free component of the forcing terms. 


3.2. Naturally Forced Receptivity. In this first subsection concerning the divergence free component 
of the forcing, the mathematically easier case of naturally forced receptivity is examined. The receptivity is 
not the result of a changing mean flow which will be studied numerically later as a homogeneous initial value 
problem, but rather the examination of a model of the forcing produced by the nonlinear interactions of 
two continuous disturbances ( e.g acoustic disturbance or free stream turbulence with surface roughness or 
surface blowing and suction). As mentioned in the introduction, this problem has two different mathematical 
formulations depending on whether the problem is viewed as a perturbation series, or the disturbances caused 
by surface effects are included in the mean flow. The main feature of this type of forcing as opposed to forced 
receptivity, examined in the next subsection, is that the forcing terms are continuous and satisfy the same 
boundary conditions as the flow at the plate, i.e the forcing and its y derivative are zero for y = 0. This 
means that the flow has evolved in such a way that there should not be any extra vorticity generated at the 
plate in order to compensate for the prescribed forcing. It is determined that in order to have this feature, 
the component B of the force must be prescribed and not A B. In the freestream and intermediate shear 
layer, the governing equations are 


(3.11) 



At; = A B, 


and 

w — sin (pU’v — A 0 . 

A Green’s function for the solution for a single harmonic in time is developed by prescribing a delta 
function source. Since the solution differs depending on whether the force is in the freestream or the 
intermediate shear layer, both are developed at once by using 

(3.13) B - r i (a, ff, fti , j/i )e iQl ( <5(y - yi) + r 2 (a, /3,U 2 ,y2)e ,n ' 21 6(y - j /•>), 


(3.12) 


~ iaU 
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with j/i > 1 so that a forcing in the freest ream is represented, and 0 < y 2 < 1 so that a forcing in the 
intermediate shear layer is represented. •> are the nondimensional frequencies of the forcing and T* •> is a 
measure of the forcing amplitudes. 

The solution is just simply 

(3.14) v - Ti «/>](<)<% - yi) + T-zi'iitWy - 1/2), 

where 


(3.15) V’i (<) = 

and 


(3.16) Mt) = | 

Since the solution (3.14) automatically satisfies all boundary conditions, there is no need for an analysis of 
the viscous sublayer at this order. 

Resonance is found to be possible when the nondimensional frequency U of the forcing is equal to the 
nondimensional wavenumber a times the value of the local nondimensional velocity. In dimensional terms, 
the resonance occurs whenever the forcing in the vorticity equations is advected with the flow. 

After integrating Squire’s equation (3.12) the general result becomes 

(3.17) w = w 0 e iaUl + sin 4> U’ f v( V ,i)e iaUit ~ T] dt - [ 

Jo Jo ft 

valid both in the free stream and in the intermediate layer. In the freestream, where U = 1 and U f = 0, 
w is independent of v. The second integral of (3.17) shows resonance in the freestream occurs provided the 
forcing time dependence varies as 

(3.18) A, C ~ e iaUt = e iat , 

the same condition shown for resonance to occur for the v velocity component. Both the magnitudes of v 
and of w increase linearly in time for this resonant mechanism. If forcing is in the intermediate layer and 
the nondimensional frequency of the streamwise and spanwise forcing is proportional to oL r (?/•>). 

(3.19) A y C - e jorJ/2< , 

then the magnitudes again increase linearly in time. However, if the frequency of the spanwise vorticity 
forcing is — otyi, the magnitude of v grows linearly in time, and the first integral of (3.17) shows the 
magnitude of w increases quadratically in time at the location of the forcing, i.e., at y — y 2 . Thus, the 
energy will increase as rapidly as t A rather than t 2 . Since any solution with fast growth might contribute 
to nonlinear interactions before the Tollmien-Schlichting waves become unstable, this resonance condition is 
important in the hunt for bypass mechanisms. 


if 

te tat if i = a, 


i« 2 f ioy 2 f 

—rrjz r— if Uo ^ 0?/‘> 

*( — aj/2) J “ 7 ~ 

te iay2i if il) = ay 2 . 
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3.3. Forced Receptivity. Again, only the divergence free component of forcing is considered and the 
governing equation is (3.11). However, instead of specifying B , one must specify A B. The consequences of 
this choice is that the analysis is much more involved since it must allow for generation of vorticity at the 
plate. Start by choosing 

(3.20) AB = ri(a,./M},)e i$2l '<5(2/ -y,) + T- 2 (a, a,U 2 )e iu ^S(y - t/ 2 ), 


with definitions of the constants the same as the previous subsection. The solution that meets the y — > oo 
boundary condition is 

(3.21) v = - H(0)e-^ u ~ l) ) e iat + ff(f)fi _ ** _1> - 

for y > 1 and 

(3.22) v = v R {y,t ) + D(t)e^ y ~ x) + E{t)e-^ y ~ X) + y - y21 , 

27 

for 0 < y < 1. The functions 0 1,2 are the same as before. Here, i?ft(y,f) satisfies the equation 

(3.23) Avr = e iay, A w 2 , f; fi (j/, 0) = ^(y) - Z)(0)e Jf< »- 1) - E(0)e-^ v ~ l) . 


The quantities iq and tb are the initial conditions for i? in the freestream and in the intermediate shear 
layer, respectively. The functions H , D and E are the coefficients of the homogeneous solution of A# = 0 
(ir rotational part of r), and must be determined by matching with the solution of the intermediate main 
shear layer with the solution in the freestream at the interface y = 1 and with the viscous sublayer at y — 0. 
Since* the time dependence is the same as before, the previous remarks about resonance and bypass also 
apply here. 

The forcing has been set in the freestream or outer reaches of the boundary layer where viscosity is 
neglected. To complement this, a uniformly valid solution that does satisfy the boundary conditions where 
v = r t/ = w = 0 at y = 0 must be established and this requires viscosity near the plate. This was not necessary 
for the case of naturally forced receptivity where the solutions in the outer and inner layers automatically 
satisfied the boundary conditions. For the case of forced receptivity, solutions from a viscous sublayer must 
be found that satisfy the boundary conditions and then asymptotically matched to the solutions of the 
intermediate layer. 

Near the plate, the governing equation for v is given by 


(3.24) 


( — iotU 1 At? — eAAt;, 

\dt ) 


where f = R~ { . 

The solution v is found by rescaling the vertical component y in equations (3.24) as Y — yj^ft (7 = 0(1) 
at most). The leading order equation formed in the limit e — ► 0 is 


(3.25) 


d ( d 2 v \ _ d A v 

dt V^rv “ W 1 ' 


A similarity solution of (3.25) for v = \Jevy , where the amplitude factor is introduced in anticipation of 
matching, is sought which leads to the introduction of 


(3.26) 


ff 2 Vy 

~0Y^ 




r) 
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A d$ n d<S> & 2 <f> n 

4T m - 2 "a^ - W = 0 


with rj = Y/2y/i as the similarity variable and r = t. Then equation (3.25) becomes 

(3.27) 

Integrating (3.26) gives 

pY pY' pi) pi) 

(3.28) t)v= / / *(Y”l2\Tt,t)dY"dY 9 = / / 4r$(r/ , ,r)df/"dr/ 

./() j 0 7o ./0 


Upon using the Mellin transform of the function <^(r/, r), defined as 

(3.29) = f r s ~ l $(7},r) dr, 

Jo 


the Mellin transform of vy and are determined to be 


or 


(3.30) 
and 

(3.31) 




44> M (t/\ s -f 1 ) dif 1 dr}\ 




= t" 2* M (ir,R+\)&i. 


At this stage, the function 4> (and therefore 4> M ) is undetermined; however, the analysis can proceed by 
noting that the plate conditions have been met and 4> A/ satisfies 


(3.32) 


$ A/ " + 2 t?* a// + 4 s4> m = 0, 


where primes denote derivatives with respect to 7 ) holding s constant. Since the solution must eventually 
match to the intermediate shear layer, the independent solution that behaves as 0 A/ ~ rf 2s must be discarded 
and 


(3.33) 


*"(»?,*) = C] (a)*, (»m) 


with $1 being the independent solution of (3.32) with <f) M ~ e. v ~ . With this solution, the quantity (^-) A/ 
is a function of s as 7/ — ► 00 . Defining 

(3.34) vr(s) 


= J~ 2C iO+ 1 )«. (*/",« + \)H 


as the Mellin transform of the unknown function V\(t), the asymptotic form as 5 — > oc of the viscous 
sublayer solution is 


(3.35) 


Vy(Y,t) ~ V’i(t)l' + Vo(f). 


The function Vo is related to V\ which is, in turn, determined by matching (3.30) with the intermediate layer 
solution (3.22). To achieve this matching requires 


(3.36) 
and 

(3.37) 


0 = v R ( 0, t) + D(t)e~i + E(t)e 7 


£2 

27 


■ 0 2 (f)c 


r 2 


-7 y-2 


V'l (t) - v Hy (0,t) + 7 D(t)e 7 - 7 E^e 1 - -rp-iiOe yy2 . 
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In like fashion the w equation can be examined in the sublayer. The governing equation is 


(3.38) 


' 8_ 
dt 


iaU 


u! — eAw = sin 4> U'v 


since there are no inhomogeneous terms due to forcing in this part of the boundary layer. By rescaling 
y = y/y/f as before and letting w = y/fW, the equation reduces in leading order to 


(3.39) 


dW d 2 W 
~dT dV 2 


sin 4> My (T, t ) 


which must be solved subject to H'(V,0) = 0 and W(0,t) = 0. Using the similarity variables and the Mellin 
transform \V M equation (3.39) becomes 


(3.40) 


W M " + 2t)W Mi + 4 sW M = 4 sin 4> ijy 1 (rj, s + 1). 


The proper matching conditions are inferred from (3.40). Specifically, as T increases, the homogeneous 
solution that grows as r/ 2s must be discarded and the homogeneous solution that decays as e~'> is kept, but 
does not affect matching. The particular solution leads to linear growth in Y dependent on the functions 
V,(0 and l'o(f). 

Assuming that the solution is started from zero initial conditions (i.e., tq = v-> = 0) leads to 
(3.41) v n (y,t) = -£)(0)e ,,!/_1) - E(0)e^«- l) . 


Combining (3.36) and (3.37) leads to 

(3.42) 
and 

(3.43) 

continuity of v results in 

(3.44) 


E(t)-E( o) = -MLe-\ 


D(t) - D{ 0) = «)^ ( f) ; 

27 27 


H(t) = H(0)e iat + + ^e- 5(1 - y2 V’2(<); 

27 ' 27 


and continuity of fi at y = 1 results in 


7 (H(t) - iaH(t)) - 

= y(^D (t) - ia( 1 - h(D(t) - T»(0))^ - 7 (e (t) - ia( 1 + b(E(t) - £(())) 


(3.45) - ye _1(1_!,2) (fa (t) - ia( 1 + ■ 

These last four equations constitute a system for the unknown functions D(t), E(t), H(t ) and Vi(<). Thus, 
the unknown function from the wall layer, Vj(f) is related to the forcing in th«> freestream and within the 
intermediate shear layer. At this point, a consistent solution for m has been constructed. The solution for 
w can also be constructed from (3.17) and automatically matches with the viscous sublayer solution. We 
note that the solutions exhibit resonance as in the case of the natural receptivity problem with the main 
difference between the two being that solutions for the forced receptivity problem immediately penetrate 
all the way to the plate even if the forcing is in the free stream. This highlights the elliptic components of 
the linear parallel disturbance equations. With this penetration into the boundary layer, it is expected that 
the forced problem leads to a much higher degree of receptivity than would the natural problem. This is 
confirmed by numerical calculations. 
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4. Receptivity and Resonance. From the previous section, we have seen that it, is possible to excite 
the continuum by a resonance like behavior. Now the question to be asked is whether or not this can lead 
to receptivity? From the work of Salwen and Grosch [10], it is obvious that any arbitrary disturbance in the 
freestream has (mathematically) a small non-zero contribution from the discrete spectrum (sE.e., Tollmien- 
Schlichting waves). This was demonstrated numerically by Lasseigne et a/., [8] where the solution to the 
initial value problem with a disturbance localized about a point yo in the freestream produced an unstable 
Tollmien-Schlichting wave. An examination of Hill’s [7] work shows that, if there is no initial disturbance, 
but rather a forcing that is turned on and later turned off, the solution of the forced problem can be re- 
expanded in terms of the linear eigenfunctions and the growth in the Tollmien-Schlichting waves calculated. 
Thus, even if the forcing is to lie strictly outside of the boundary layer, it must still excite an unstable 
Tollmien-Schlichting wave since there is little possibility that the resulting solution is orthogonal to the 
adjoint of the Tollmien-Schlichting wave. This is an example of forced receptivity that does not rely on 
any disturbance within the boundary layer proper or on any degree of non-parallelism in order to excite 
the unstable Tollmien-Schlichting wave. The numerical calculations presented in this paper characterize the 
strength of this mechanism. 

That this receptivity mechanism does not require non-parallelism does not mean that. 11011 -parallelism 
can be neglected when considering receptivity. After the forcing has been turned off at a finite time, it 
is not physically correct to expand in terms of the initial eigenfunctions since the disturbance resulting 
from the forcing should be considered as the initial conditions at a higher Reynolds number. Thus, even 
the continuum modes that have been excited by the forcing term at the original Reynolds number could 
produce a Tollmien-Schlichting w r ave at the higher Reynolds number. However, the converse may also be 
true. Mathematically, this is a result of the eigenfunctions at one Reynolds number not being orthogonal 
to the adjoint eigenfunctions at another Reynolds number. Numerical calculations presented in this paper 
wall show T that through this mechanism, non-parallel effects downstream of the leading edge and even into the 
unstable region play a major role in the receptivity problem . In particular, Nishioka and Morkovin [9] note 
that for the forced receptivity problem Tollmien-Schlichting wave packets demonstrate growth far exceeding 
the growth of a theoretical wave packet in which it is assumed that no further seeding of the Tollmien- 
Schlichting waves takes place upstream of branch I of the neutral curve. Others have found such phenomena 
when considering naturally forced receptivity to freestream disturbances. It is also noted here that growth 
rates matching the theoretical values are found for the vibrating ribbon problem in which the continuum 
plays much less of a role. 

The development of the theory used here for both parallel and non-parallel boundary layers is straight- 
forward and relies on the assumption that the parallel theory is at least locally applicable for each value of 
the Reynolds number, i.e., the same assumption used to derive the Orr-Sommerfeld equation and as a basis 
for analyzing the effects of non-parallelism on a single Tollmien-Schlichting w v ave. Assume that the governing 
equation (2.13) for the parallel temporal problem has been written as 

(4-1) L r (v) = /(y, 


w T ith f(y, 0) = 0 and f(y,t) = 0 for t > T, The frequency cj/ is considered as the primary frequency of the 
forcing. 

Since the eigenfunctions of the Orr-Sommerfeld equation form a complete set, the solution at t = T 
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found bv solving (4.1) subject to zero initial conditions and Reynolds number R = R] is 


(4.2) 


N POC 

(y, T) = V' aj<t>i(y, R \ ) + / A(k)<t>(y r fc, ) dk 

i= i Ju 


where it is chosen so that the discrete eigenfunctions are normalized to have unit energy and the continuum 
are also normalized with respect to the energy, i.e., 


(4.3) 


^ = i 

7“ Jo 


and 

-j rk-\-d rOC 

(4.4) -77U / / / <F(v,k",Ri)(0 2 -i?)<Kv,V,Ri)dydk'dk" = l. 

4tf“7~ J k ._ s J k~s Jo 

If we choose $ as the properly normalized eigenfunctions of the adjoint solution, i.e., 

(4.5) I" - f)<Pj{y,R\)dy = 

7 Jo 


and 

(4.6) r **(y,*", JliHfl 2 - ^WyA-'MOdy = 6(k" - k% 

7“ Jo 

then the coefficients in (4.2) are 

(4.7) «, r ^(y^0(D- - J 2 )v p (y, T) dy 

7 Jo 

and 

(4.8) Mk) = -^f $*(y,k,Ri)(D' 2 - T)v p (y,T) dy. 

7“ Jo 

The coefficients a, are the increase in the amplitude of the Tollmien-Schlichting waves due to the forcing 
and .4(A;) is the increase in the continuum modes due to forcing. Since there is no forcing for t > T, the 
evolution of the solution beyond t = T is given by 


(4.9) v(y,t) = TateW'W-T'MV'Ri) + H A(k)e^e-^+ k ' 2 ^- T U(y,k,ROdk 

7^i Jo 

Equation (4.9) represents the exact solution for a purely parallel flow subjected to forcing. However, the 
boundary layer thickens as time progresses, and for a more realistic solution the effects of this must somehow 
be included. This thickening of the boundary layer is of course a continuous process, but insight can be 
gained by re-expanding the solution at t = T in terms of the eigenfunctions for a larger Reynolds number. 
For this calculation, great care must be taken with the normalizations. First, all scalings in the {/-direction 
remain as the standard scalings with R = R \ . The eigenfunctions and adjoint eigenfunctions for R= R- 2 are 
generated by replacing U(y), U'{y), and U"(y) in equations (2.13) and (2.14) with 

(4.10) U(y) -4 U(y/y/^), U'(y) -4 -±=U'(yU^), and U"(y) -4 —U"(y/^), 

where ^/xq — R 2 /R 1 • The same normalization and orthogonalization relations are assumed to apply to the 
new set of eigenfunctions. Expanding in terms of the eigenfunctions <j>(y, J?, 2 ) and <j){y, k, R?) instead of the 
expansion (4.2) gives 


(4.11) 


v{y,T) ~ R 2 ) + [ B(k)4>(y,k,R2)dh 

i= 1 
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with 


(4.12) 


>, = -/ 


^(y,R 2 )(D 2 



(y,Ri)+ [ A(k)<p(y,k,Ri)dk \ dy. 
Jo 


Taking the extreme example of a* = 0, ie., the forcing producing no Tollmien-Schlichting waves at R = i?i, 
the coefficient bj 0 since $*(?/, i?>) is not orthogonal to <j){y,k, R\) and Oj{y, R\ ) when j ^ i. This is 
a receptivity mechanism due to non-parallelism alone since no further disturbance within or without the 
boundary layer is required to initiate the gain in amplitude of the Tollmien-Schlichting wave. At first, it 
might be tempting to discount this route to receptivity. Upon examination on the inter-relationship between 
(pj(y,Ri) and $*(?/, # 2 ) one finds 

(4.13) -4 r <t>'i{y,R 2 ){D 2 -7 2 )0i(»,«i)rf»« 1 

7“ Jo 

for a wide range of R] and It is therefore natural to assume that the orthogonality relations are also 
approximately preserved and that 6,- would automatically be quite small if a* — 0. To determine whether the 
orthogonality relations are approximately preserved with changing Reynolds number, the entire spectrum 
of the Orr-Sommerfeld equation and its adjoint must be calculated along with the integrals indicated in 
(4.12). For the continuum, this is not a feasible undertaking. However, by numerically solving the initial 
value problem, this mechanism is quantifiable. The results of such a calculation are given in Section 5.2. 
Clearly, the analysis presented above is possible since the eigenfunctions of the Orr-Sommerfeld equation 
form a complete set and any function of y satisfying the same boundary conditions can be expanded using 
this set. This leads to a more direct solution of the forced problem. Considering the component of forcing 
F = (.4,Z?,C) with V • F — 0, then equation (2.13) becomes 


(4.14) 



At? + iaU”v = R~ l A At? + A B, 


with B(y,t) being the Fourier transform of the normal component of the force F that is assumed to satisfy 
B(Q,t) — B y (0^t) — 0. If it is also assumed that v = 0 and B = 0 for t < 0, then the Fourier transform in 
time to (4.14) gives 


(4.15) 


[iu — iaU] A?; + iaU n v = R 1 A An + AB 


where 

1 C°° 

(4.16) v(a,y,p,Q) = — / v(a,y,0,t)e luJt dt. 

j — 00 

Both of the time transforms v and B can be expanded in terms of the complete set of functions as 


(4-17) 


v(y) 


N no 

= ai<j>i(y) 4 - / 

f=l 


A(k)<f>(y , k) dk 


(4.18) 


B{y,Q) = y'p i {u)<f>j(y) + / ,k)<j>(y,k)dk. 

,=i Jo 


Since the eigenfunctions 0j(t/), i = 1 ... N and (j){y , fc), k = (0, 00 ) satisfy 
(4.19) iaUA<j>* — iaU”<fi * + R~ l AA<j)* = icu+Afi*, 
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there is a one-to-one correspondence between the coefficients of the expansion for v and B. Thus, the exact 
solution is written as 


(4.20) 


N 

v{y,ti) = 

1 = 1 


i(w — 



0(k,v) 

i(uj — uj(k)) 


cj)(y, k) dk. 


Then, inverting the Fourier transform in time by integrating along a contour that lies below all poles 
provides the coefficient of the Tollmien-Schlichting wave produced by the prescribed forcing, i.e.. 


(4.21) ft; = 2nf3j(u>i). 

Although the temporal aspects of this result are not necessarily new (see Hill [7], Nishioka and Morkovin 
[9], and Crouch [4] for spatial equivalents), the dependence of this receptivity factor on the vertical profile 
of the divergence free component of the inhomogeneous terms is new. Some clarification of the temporal 
results implied by (4.21) is necessary to put the calculations that follow into context. First, the non-zero 
Tollmien-Schlichting wave generated by forcing is not the result of a true mathematical resonance but instead 
is the normal response of the solution of the homogeneous linear system to an initial value problem with an 
inhomogeneity. If forcing were to be of a single real frequency, u Jj say, started at t — oc instead of t ~ 0, then 
is a delta function, and the evaluation at the complex frequency ujj of the ith Tollmien-Schlichting wave 
is identically zero. However, causality is crucial to the calculations that follow. For single real frequency 
forcing starting at t = 0, ft(u/) has a simple pole at u) f . It is the smallness of the growth rate, i.e., the 
imaginary part of that produces a near resonance like behavior for the initial value problem if the real 
frequency jJj matches the real part of u>i in the same manner as the near resonant behavior in the weakly 
damped oscillator. However, prescribing a step function forcing is not a valid model for receptivity. The 
step function is known to excite all frequencies (and hence a pole rather than a delta function behavior of 
u;)), and it is no surprise that a Tollmien-Schlichting wave should be excited. 

There is, however, a greater chance of resonance like behavior for the continuum modes of the Orr- 
Sommerfeld equation. First, it is recognized that /*; (u;) has a pole at the real frequency u ;/ for any forcing 
with pure sinusoidal behavior for t > t 0 and to » 1, no matter how smoothly the forcing changes from a 
zero value at t = 0 to this behavior at t = t. 0 . This pole has a factor of which would be small for 

the unstable Tollmien-Schlichting waves. However, the continuum inodes with u, i(k) — a + if (a 2 + k 2 ) has 
an imaginary part very near zero. For R = 1000 and a = .25, the imaginary part of cj( 0) has a value of 
.0000625 while the imaginary part of the unstable Tollmien-Schlichting wave has a value of .00302. Thus, if 
the forcing is advected with the flow, i.e., ujf = a, a very near resonance with the continuum should appear, 
just as predicted by the model problem in Section 4. 

Through (4.21) it is seen that the forced problem can be directly related to an unforced initial value 
problem. Assuming that the forcing is written as the product B = f(u — uif)G(y), then the strength of the 
Tollmien-Schlichting wave is proportional to 

(4.22) -- / *;(y)(D 2 -j 2 )G(y)dy, 

7“ Jo 

which is the coefficient of the Tollmien-Schlichting wave for the initial value problem with v(j/,0) = G{y). 
In both Salwen and Grosch and in Hill, examples are explored by replacing the Laplacian of G with a delta 
function centered at y = y 0 - In Salwen and Grosch this is equivalent to specifying the initial vorticity and in 
Hill this is equivalent to specifying the vorticity source. Since 4> * ( i; ) behaves as e _ ™ as // -4 oc, one expects 
that the the response to a localized source in the freestream would be proportional to e~~' yo . However, 


16 



the results of the numerical calculations presented later do not necessarily have this property. Instead the 
response sometimes decreases at a rate much faster than r“™°. The discrepancy is resolved by realizing 
that the specification of arbitrary initial vorticitv or an arbitrary vorticity source as a delta function ignores 
where these quantities come from and the special nature of the integral (4.22). One interpretation of the 
initial value problem is that the initial values are assumed to come from the solution of the flow equations 
at an earlier or upstream position (naturally forced receptivity). This concept is pushed upstream until 
it is no longer possible to ignore the leading edge effects. For this interpretation, the initial vorticity has 
to be consistent with a realizable velocity field and the delta function is not. The easiest way to resolve 
this is to specify the initial velocity and then calculate the initial vorticity from this specification. If this 
is done, integration by parts of (4.22) shows that there is not an e~^ yo behavior of the strength of the 
generated Tollmien-Schlichting wave for this mode of receptivity. This is as predicted by the model problem 
in Section 4. 

5. Numerical Solutions of the Linear System. For the analysis that follows, F is chosen to have a 
localized structure in the y variable so that the effects of forcing the boundary layer at various locations are 
determined. Since the problem is linear, the response to these localized forces can be considered as a form 
of a Green's function for the problem. 

The partial differential equations (2.13) and (2.14), together with the boundary conditions (2.16) at the 
plate and boundedness in the freestream, were solved numerically by the method of lines. This is a convenient 
numerical method that has has worked well in the past; other techniques are possible. The spatial derivatives 
were center differenced on a uniform grid and the resulting system integrated in time. 

In using such an approach on an unstable system of equations in order to quantify the instability, there 
is always the possibility that the instability observed in the numerical solution is actually the buildup of 
numerical errors rather than the instability resulting from the initial conditions or the imposed forcing. This 
work is particularly subject to such a criticism since it is proposed that the forcing lie entirely outside of the 
boundary layer. However, in previous work (Lasseigne et.al . , 1998), numerical errors within the boundary 
layer did not effect the solution to the initial value problem when disturbances originating outside of the 
boundary layer were considered. As the center of the initial disturbance moved away from the boundary 
layer, tighter tolerances were necessary to keep the numerical errors that accumulate inside the boundary 
layer within the proper bounds. To circumvent the issue of numerical errors producing the instabilities 
observed in the receptivity calculations presented below, all solutions presented are subject to the initial 
condition 

,c JN - y 2 e~ y ' /2 

/o°° (& 4 + (!/a 2 ) (j/ 2 (2 - y 2 ) 2 )) e~v 2 dy 

where the normalization is chosen so that the energy is unity at the initial point. This initial condition 
produces the quantifiable, unstable Tollmien-Schlichting wave seen in Figure (5.1) labeled by Bq = 0. To 
this initial value problem, forcing is added in the governing equations. It is always chosen that the forcing 
is identically zero at t = 0 and increases smoothly as a function of time. Receptivity results are presented 
as a gain in the amplitude to this unstable Tollmien-Schlichting wave. The drawback to including an initial 
condition in the forced problem is that by having two disturbances in the solution, the time marching of 
the numerical scheme must be slower to resolve the entire solution. However, efficiency of the numerical 
scheme is not of issue for this work, but proper modeling is. Even at its slowest, numerical integration of 
the linear, temporal stability equations is considerably faster than direct numerical simulation of the spatial 
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PlG. 5 . 1 . Response to slowly ramped single frequency forcing applied at yo ~ 6 . Frequencies are u/ = 0 . 05 , 0 . 25 ,. Dashed 
line for u — .25 with foreing stopped at t — 1000. 


problem. Thus, this remains an effective tool for exploring the nature of the solutions under a wide variety 
of conditions and guiding the more expensive numerical simulations. 


5.1. Parallel Calculations. The first issue addressed by numerical investigation is whether or not the 
resonant behavior, as described in previous sections, actually occurs for the Blasius flow. To this end, the 
forcing function B is chosen and A and C are such that there is no forcing in the pressure equation nor 
in the normal vorticity equation. A general purpose function that satisfies the boundary conditions and is 
proportional to a delta function centered at yo in the limit as a y — > 0 is 


(5.2) 


B(y) = B 0 r(t)e 


(y/y o) 2 e 


2f>-(y-yof/cr~ 


{S^liylyoYc-^^/^Ydy) 


1/2 


The Laplacian of this function is used in the equation (2.13). 

The time histories of the perturbation energy using the parameters values B 0 = 1000, R = 1000, 
a = 7 = .25, 0 = 0, o y = .5, and y 0 = 6 are shown in Figure (5.1). The forcing frequencies are ujj = 0.05 
and Uf — 0.25 (for comparison, the Tollmien-Schlichting frequency is .0874). The function r(£) = 1 — t a * , 
o t = 100, provides a smooth, slow increase in the forcing. After a short transient, the constant energy curve 
when u>f ^ a represents the single frequency particular solution of the forced governing equations. The 
energy remains constant until the energy of the forced solution and the energy of the unstable Tollmien- 
Schlichting wave introduced by the initial conditions are of the same size. The curve with u)f = a shows 
greater growth in energy than the curve with ujf ± a, and this curve does not level off to a finite value. The 


18 





Fig. 5.2. Response to transient forcing applied at yo — 1.0, 3.0, 5.0. Frequencies are uij = 0.05, 0.25. af 


10, 25, 50. 


growth in energy for this resonant frequency approaches t 2 but is slightly less owing to damping by viscosity. 
This is as predicted in the inviscid analytical solutions. 

In these calculations, no formal receptivity is expected since the imposed conditions are chosen to 
closely approximate single frequency forcing. However, according to the theory presented in Hill [7] and in 
Section 4, if the forcing is removed, then an unstable Tollmien-Schlichting wave must be generated. The only 
remaining question is the size of this wave. The dashed line in Figure (5.1) shows the results of smoothly 
removing the forcing at t = 1000. A significant Tollmien-Schlichting wave is not generated by forcing 
external to the boundary layer at this resonant frequency even though the energy of the disturbance shows 
tremendous growth prior to removing the forcing. Mathematically, the Tollmien-Schlichting wave must have 
been generated, but, since the conditions chosen (yo large and r(t) smooth ) make (ii(uj) an extremely small 
number, it does not appear in Figure (5.1) . 

Transient effects of forcing are explored by calculating the response to the forcing with 


(5.3) 


r(t) = 


/ , \ 2 

f L I e -(t-t 0 )' 2 /a; 

VoJ 


Figure (5.2) shows the time histories of the perturbation energy for eighteen cases: to = 200, Bo — 10 6 , 
a = 7 = .25, <p = 0, (T y = .5, erf = 10, 25, 50, u;/ = 0.05, 0.25 and yo = 1.0, 3.0, 5.0. Even though the forcing 
is applied well outside of the boundary layer (up to yo = 5), there is an increase in the Tollmien-Schlichting 
wave amplitude owing strictly to this forcing. Clearly, the dominant parameter is the vertical position of 
the localized forcing function. This might even be unexpected since the Tollmien-Schlichting wave and its 
adjoint both are proportional to exp (— 71 / 0 ) at the point of forcing. This factor cannot account for the more 


19 




that five orders of magnitude decrease seen between groups of solutions as the parameter y 0 varies. Also of 
important note in Figure (5.2) is the slow algebraic decay after cessation of the forcing when y {) = 5. This 
is the hallmark of the response as predicted in Grosch and Salwen [10] for the continuum modes. 


It is seen above that even for a single wavenumber 7 = .25 and a single Reynolds number R — 1000, 
the strength of the receptivity response depends on a wide number of parameters forcing introduces into the 
problem. Using a forcing function of the form (5.2) subject to (5.3), the free parameters for this one point on 
the Orr-Sommerfeld diagram are: # 0 , fen uj/, u 1 1 and y$. Bo is just a scaling factor in a linear problem, 

t 0 should be chosen sufficiently large such that forcing does not start immediately, and there is relatively 
weak dependence on ujf and a t (as compared to ;t/ 0 ). It is interesting that these two parameters do not have 
more of an affect considering that the width of the forcing spectrum is solely determined bv them. 

From the theory in Section 5, the receptivity factor is v(u>'ts ~~ ^f) (where f is the Fourier transform 
of (5.3)) times the coefficient of an eigenfunction expansion of the function B(y) can be measured by ex- 
trapolating the time histories of E(t) back to t = 0 . If a T $ is the complex coefficient of the unstable 
Tollmien-Schlichting wave normalized as in (4.3). then the intercept of the straight line extrapolation of the 
Tollmien-Schlichting energy growth seen in Figure (5.2) represents |«ts| 2 - Thus, the gain in magnitude of 
the quantity 


(5-4) 


Ats = 


*TS 


T a 


TS\ 


l TS \ 


\Bo\ 


will be referred to as the receptivity factor and used to characterize the response to the forced problem. In 
the limit |J 3 0 | 00 , the contribution to the receptivity factor is from the forcing only and not from the 

initial condition. This coefficient has been determined for the eighteen cases included in Figure (5.2) and is 
tabulated in Table 1 with the Orr-Sommerfeld coefficient determined by integrating the disturbance equations 
with v(y, 0) = B(y) and no forcing. Also included in Table 1 is the value of r(u ;ts — &f) which has been 
determined numerically. In every case, when r{uTS ~~ u; /) I s multiplied by the coefficient determined by the 
initial value problem, the theory agrees with the forced calculations. Additional values of v{u)rs - Re(u ^ ts))> 
i.e. forcing at exactly the frequency of the Tollmien Schlichting wave and the predicted receptivity factor 
are calculated and tabulated with the previous results. For these parameter values, the tuning of the forcing 
to match the Tollmien-Schlichting wave does not have much affect. 

In the previous calculations, the function B has been specified and therefore these curves represent the 
case of naturally forced receptivity. In Figure (5.3), the same function was used to specify A 2 B and the 
differences in Ats as a function of y 0 is shown. As the disturbance location moves toward the freestream, the 
immediate generation of Tollmien-Schlichting waves have strengths proportional to e~ lv as was previously 
predicted (see Hill, 1995). The strength of immediate generation of Tollmien-Schlichting waves for the 
naturally forced receptivity drops off at a much faster rate as the forcing location moves toward the freestream. 
Thus, receptivity to turbulent disturbances which are shielded from the plate region by mean flow shear will 
have less than expected strength considering that the Tollmien-Schlichting waves decay as e ~ ly . 
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Table 1 

Receptivity factor for transient forcing. 


*) 

Uf 

I/O 

r(cjrs - uj) 

Initial Value 

Forced 

10 

.05 

1 

3.052 

1.195* 10* 

3.647* 10 1 

10 

.05 

3 

3.052 

3.664 * 10- 2 

1.118* 10" 1 

10 

.05 

5 

3.052 

1.162 * 10~ 5 

3.539 * 10~ 5 

25 

.05 

1 

4.800 

1.195* 10 1 

5.736* 10 1 

25 

.05 

3 

4.800 

3.664 * 1()-' 2 

1.759* 10-’ 

25 

.05 

5 

4.800 

1.162* nr 5 

5.567 * 10“ ft 

50 

.05 

1 

6.730 

1.195* 10' 

8.041 *10' 

50 

.05 

3 

G.730 

3.664 * 10 2 

2.466* 10- 1 

50 

.05 

5 

6.730 

1.162* 10" 5 

7.804 * 10- 5 

10 

.25 

1 

2.867 

1.195* 10 1 

3.425* 10 1 

10 

.25 

3 

2.867 

3.664 * 10~ 2 

1.050* 10“ 1 

10 

.25 

5 

2.867 

1.162 * 10— 5 

3.324 * 10- 5 

25 

.25 

1 

4.106 

1.195 * 10 1 

4.906* 10 1 

25 

.25 

i 

3 

4.106 

3.664 * 10“ 2 

1.504* 10-' 

25 

.25 

5 

4.106 

1.162* 10 5 

4.761 * 10- r> 

50 

.25 

1 

4.924 

1.195* 10' 

5.883 * 10 1 

50 

.25 

3 

4.924 

3.664 * 10“ 2 

1.804 * 10 _1 

50 

.25 

5 

4.924 

1.162* 10“ 5 

5.710 * 10- 5 

10 

.08744 

1 

3.063 

1.195* 10 1 

3.660 * 10' 

10 

.08744 

3 

3.063 

3.664 * nr 2 

1.122* 10" 1 

10 

.08744 

5 

3.063 

1.162* 10“ 5 

3.559 * 10- 5 

25 

.08744 

1 

4.843 

1.195 *10* 

5.787* 10 1 

25 

.08744 

3 

4.843 

3.664* 10~ 2 

1.774* 10- 1 

25 

.08744 

5 

4.843 

1.162* 10 5 

5.628 * 10“ 5 

50 

.08744 

1 

6.849 

1.195* 10 1 

8.184 * 10 1 

50 

.08744 

3 

6.849 

3.664 * nr 2 

2.509* 10“* 

50 

.08744 

5 

6.849 

1.162* 10” 5 

7.958 * 10- H 


5.2. Non-parallel Calculations. In the standard Orr-Sominerfeld approach, the value of R corre- 
sponds to a specific downstream position and therefore to a particular thickness of the boundary layer. The 
Orr-Sommerfeld solution approach also focuses on a single mode (the least stable if R < R c or the unstable 1 
mode if R > R c ) which is assumed to be periodic in the streamwise direction. The basic non-parallel theory, 
as set out in Gaster [5], is properly considered as a problem of a spatially evolving flow. However, it appears 
that what is considered by this theory is the evolution of only one mode and not the interaction of other 
modes. Thus, it is still unclear how an individual disturbance, written as a sum of all eigenfunctions, evolves 
as the value of R changes smoothly. The techniques employed here allow for the unique opportunity to 
measure the linear interactions between the eigenfunctions at two different downstream positions. 

The response to forcing of the form (5.2) is shown in Section 5 to be linked to the solution of the initial 
value problem and this result is confirmed by the numerical experiments for the parallel problem. Therefore, 
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>0 


Fk;. 5.3. Difference in receptivity factor between forced modes (diamonds) and naturally -forced modes ( squares ) of recep- 
tivity as a function of disturbance location yo- 


for the purpose of clarity, non-parallel effects are explored through a series of initial value problems, which 
are directly applicable to the problem of natural receptivity; the implications for the naturally forced problem 
can be inferred. In the spirit of using the temporal problem to explore a problem that is properly considered 
spatial, the Reynolds number in the disturbance equations (2.13-2.14) will be held constant and the mean 
How will be amended to represent different downstream positions. The transformation is given by (4.10). The 
standard Orr-Sommerfeld diagram represents holding the mean How constant (x 0 = 1) and the downstream 
evolution of a single wavelength disturbance traces a ray emanating from the origin of a R-a diagram. In this 
approach (/? is held constant and j: 0 varied), the downstream evolution of a single wavelength disturbance 
traces a horizontal path in the aro-o diagram which is shown for R = 1000 in Figure (5.4). 

In the first series of calculations, an initial value of the form (5.2) with r(f) = 1 and up = 0 is added to 
the benchmark initial disturbance (5.1) and the receptivity factor T is computed. Figure (5.5) shows Ats 
as a function of y 0 . The receptivity factor is found to be decreasing slightly faster than exponentially for 
increasing yo as previously mentioned, but it is also shown to be an increasing function of x 0 (by as much 
as two orders of magnitude when yo =4). When the disturbance is inside the boundary layer for all values 
of Xq {yo = 1), the increase in receptivity factor as x 0 increases is about half an order of magnitude. All 
values shown are for the unstable region of the Orr-Sommerfeld diagram where it is possible to calculate the 
receptivity factor using this approach. For y 0 = 5 the disturbance is completely outside of the boundary 
layer and the increase in the receptivity factor is now several orders of magnitude (but the factor is still 
small). The solutions for smaller values of x 0 when yo = 5 could not be properly resolved using the same 
parameter values, but the extrapolation of the results is obvious. 
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In the next series of initial value problems shows order of magnitude increases in the receptivity factor as 
j*o increases and y G large are due to the continuum modes and other Tollmien-Schlichting modes from lower 
value of £o feeding the Tollmien-Schlichting modes at the higher x G , For x 0 = 0.8, an unstable Tollmien- 
Schlichting mode is found by calculating the solution to the governing equations for large time. This mode 
is then normalized to have unit energy and is used for an initial condition. Figure (5.6) shows that the 
variation of the Tollmien-Schlichting amplitude as x 0 varies from stable through the unstable region of the 
Orr-Sommerfeld diagram is slight. The curve for xq = 0.4 starts out with exponential decay which is to be 
expected if the initial condition is mostly the stable Tollmien-Schlichting wave and then decays algebraicly 
as expected since only the continuum is left for large times. All curves extrapolated back to t = 0 along 
the exponential part of the growth or decay give values of receptivity factors between 0.8 and 1.25. Next, 
the solution for initial values of the form (5.2) with ~ 3 and xo — 0.8 is calculated twice: once as is and 
once with an additional Tollmien-Schlichting wave added that is amplitude and phase matched to nearly 
cancel the instability. The results are seen in Figure (5.7) where the divergence of the solutions clearly show 
that the value of ars for the second case is near zero. The solution using the same two initial values an' 
calculated with x 0 = 1 and x 0 = 1.2. Surprisingly, it is seen that a majority of the Tollmien-Schlichting wave 
at these higher values of x () do not come from the Tollmien-Schlichting wave at x 0 = 0.8 but rather from 
the part of the solution that produces no Tollmien-Schlichting wave at x 0 = 0.8. The results here are rather 
ominous. If then' is any additional disturbance in the outer edges of the boundary layer or near freest ream, 
these disturbances feed directly into the Tollmien-Schlichting wave and will produce a growth rate greater 
than (sometimes very much greater than) the predicted value, even when that predicted value accounts for 
all of the non-parallel effects associated with a single Tollmien-Schlichting wave. 

6. Conclusions. It has been shown that the techniques previously developed by the authors to in- 
vestigate various aspects of the temporal stability problem can also be applied to investigate the problem 
of receptivity. It is shown that resonance with the continuum can occur, and this must, be considered 
when investigating bypass mechanisms. Also shown is that the form in which the forcing function is intro- 
duced into the governing equations has great significance when determining the strength of the generated 
Tollmien-Schlichting wave. Perhaps most importantly, it is shown that the transfer of the solution from one 
downstream location to another is not one-to-one in terms of the eigenvalues, and that the continuum at one 
position feeds into the Tollmien-Schlichting wave at another. 
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